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Abstract. Scattered light images of the optically thin dust disk around the 5 Myr old star HD 141569 A have 
revealed its complex asymmetric structure. We show in this paper that the surface density inferred from the 
observations presents similarities with that expected from a circumprimary disk within a highly eccentric binary 
system. We assume that either the two M stars in the close vicinity of HD 141569 A are bound companions 
or at least one of them is an isolated binary companion. We discuss the resulting interaction with an initially 
axisymmetric disk. This scenario accounts for the formation of a spiral structure, a wide gap in the disk and a 
broad faint extension outside the truncation radius of the disk after 10-15 orbital periods with no need for massive 
companion(s) in the midst of the disk resolved in scattered light. The simulations match the observations and the 
star age if the perturber is on an elliptic orbit with a periastron distance of 930 AU and an eccentricity from 0.7 to 
0.9. We find that the numerical results can be reasonably well reproduced using an analytical approach proposed 
to explain the formation of a spiral structure by secular perturbation of a circumprimary disk by an external 
bound companion. We also interpret the redness of the disk in the visible reported by Clampin et al. (2003) and 
show that short-lived grains one order of magnitude smaller than the blow-out size limit are abundant in the 
disk. The most probable reason for this is that the disk sustains high collisional activity. Finally we conclude 
that additional processes are required to clear out the disk inside 150 AU and that interactions with planetary 
companions possibly coupled with the remnant gas disk are likely candidates. 
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1. Introduction 

. Asymmetries and annular structures are common obser- 
' vational features regardless of the emission process (scat- 
tered light or thermal emission) for the handful of gas- 
free and optically thin disks currently resolved around 
Main-Sequence stars. Surface brightness maxima peak far 
from the star from dozens up to about one hundred AU. 
Attempts to explain radial and azimuthal structures in- 
volve massive un- resolved planet (s), trapping dust parti- 
cles into resonances due to either radial migration of par- 
ticles sensitive to Poynting-Robertson drag and/or to ra- 
diation pressure (e.g. eEri: Ozcrnoy et al., 2000; Quillen 
& Thorndike, 2002) or to outward migration of the planet 
(e.g. Fomalhaut and Vega: Wyatt & Dent, 2002; Wyatt, 
2003). Other attempts involve external stellar compan- 
ion(s) either bound and observed (e.g. HR4796: Wyatt 
et al., 1999) or unbound and currently unobserved but 
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having recently approached the close vicinity of the disk 
(flyby scenario, e.g. /3Pictoris: Larwood & Kalas, 2001). 
In the edge-on disk of /3 Pictoris, the vertical asymmetries 
are explained by the precession of planetesimal orbits in- 
duced by an inner planetary companion on an orbit that 
is inclined to the dust disk (Mouillet et al., 1997) and the 
effect of radiation pressure acting on the smallest grains 
(Augereau et al., 2001). But note that the inner planet 
in this model can be replaced by any inner misaligned 
mass distribution with the appropriate quadrupole mo- 
ment components and is not dependent on any partic- 
ular planet postulate. An alternative scenario involving 
dusty clump formation through stochastic collisions be- 
tween large planetesimals has been proposed by Wyatt 
& Dent (2002) to explain the asymmetries noticed in the 
Fomalhaut disk. 

We explore in this paper a source of asymmetry for 
the optically thin dust disk surrounding HD 141569 A, a 
B9.5V-A0V star located at about 100 pc according to 
Hipparcos measurements. Coronagraphic images from the 
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Fig. 1. HST/STIS visible image of the optically thin dust disk 
around HD 141569 A from Mouillet et al. (2001). The two M 
companions, HD 141569 B and C, located in the North- West 
region lie outside of the image (see text for precise location). 



visible to the near-infrared have revealed the complex mor- 
phology of the dusty circumstellar environment of this old 
Herbig star (Augereau et al., 1999; Weinberger et al., 1999; 
Mouillet et al., 2001; Boccaletti et al., 2003; Clampin et 
al., 2003). One can construct the following sketch of the 
shape of the dust disk beyond 100-120 AU (the edge of 
the coronagraph) as it appears in scattered light (see also 
Figure 1) : 

— the disk is composed of two annuli peaked at 200 AU 
and 325 AU from the star with centers shifted by 
20-30 AU almost along the minor axis of the disk 
(East- West direction). The outer ring actually shows 
a tightly-wound spiral structure, 

— the two bright predominant annuli at 200 AU and 
325 AU have between them a darker ring or "gap" . 
This gap is radially wide compared to the two annuli, 

— the two bright rings show out of phase brightness 
asymmetries of up to factors of 2.5 3 for the outer ring 
in the visible. These asymmetries cannot be explained 
by invoking scattering properties of the dust grains, 

— an extended diffuse emission is present in the North- 
East of the disk and is detected up to more than 
600 AU, 

— the disk brightness sharply decreases between 200 and 
150 AU rapidly reaching the background level of scat- 
tered light images interior to 150 AU. This behavior is 
suggestive of a strong, but likely not complete, deple- 
tion of dust inside 150 AU, 

— structures with smaller spatial scales are also present 
such as a radially thin arc superimposed on the dark 
lane/gap at a distance of 250 AU from the star, 

— the maximal vertical optical thickness of the outer ring 
is ^ 2 ± 1 X 10~^ in the visible and the near-infrared. 



The dust content inside 100-120 AU remains almost 
totally unconstrained despite marginally resolved images 
in the mid- infrared (Fisher et al., 2000; Marsh et al., 
2002) which indicate a confirmation of dust depletion 
inside ~100AU (Fig. 4 from Marsh et al., 2002). A total 
midplane optical depth in the visible of ^0.1 has been 
estimated by Li & Lunine (2003) indicating that the disk 
is optically thin in all directions. 

HD 141569 A is not isolated but has two low-mass stel- 
lar companions HD 141569 B and HD 141569 C located at 
7.54" and 8.93". Their position angles (PA) are 311.3° 
and 310.0° respectively (Augereau et al., 1999). We show 
in this paper that the gravitational perturbation of the 
HD 141569 A disk by the detected stellar companions gives 
a natural explanation for some of the broadest features ob- 
served in scattered light as long as one of the companions, 
or both if bound, is on an orbit with high enough eccen- 
tricity. 

We detail in section 2 our motivations for exploring 
the impact of the observed companions on the shape of 
disk and we give a description of the dynamical model we 
used to address this issue. The numerical results shown in 
section 3 are compared with an analytic solution to the 
problem in section 4. A surface density consistent with 
the resolved images of HD 141569 A is obtained in section 
5 and we discuss the implications for the dynamics of the 
companions. In section 6, we interpret the redness of the 
disk in the visible measured by Clampin et al. (2003) in 
terms of minimal grain size in the disk and we discuss 
the consequences of these results. We finally point out the 
limitations of our dynamical approach in section 7 and 
indicate directions for future work. 

2. Dynamical modeling 

2.1. Motivations 

We discuss in this paper the possible gravitational influ- 
ence on the dust disk of the stars HD 141569 B and C ob- 
served in the close vicinity of HD 141569 A. Motivations 
for exploring the potential gravitational influence of the 
detected stellar companions are threefold: 

1. Based on astrometric and radial velocity measure- 
ments, Weinberger et al. (2000) argue that B and C, 
identified as M2 and M4 pre-Main Sequence stars re- 
spectively, have a high probability of forming with 
HD 141569 A a bound triple system. The three stars 
have similar ages : 5±3 Myr for the two M companions 
(Weinberger et al., 2000) and ~5Myr for HD 141569 A 
(Merm et al., 2003). If HD 141569 A, B and C are in- 
deed gravitationally linked, B and/or C are likely to 
perturb the disk surrounding the primary and are nat- 
ural candidates for structuring the disk. In any case the 
possible effects of these observed probable companions 
has to be explored before invoking other sources of 
asymmetries in the disk. Moreover the very latest im- 
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ages of the disk with HST/ACS strongly support our 
approach (Clampin et al., 2003). 

2. Could undetected compamon(s) in the midst of the 
disk explain the gap at 250 AU? Weinberger et al. 
(1999) show that a companion with a mass larger than 
3 times that of Jupiter would have been detected by the 
HST/NICMOS instrument but on the other hand an 
object with a mass 1.3 times that of Jupiter would be 
sufficient to form and clear the 250 AU region. While 
the presence of such companion(s) can not be ruled 
out, current planet formation models do not readily 
account for the in situ formation of Jupiter-like objects 
at a distance of a few hundred AU around a ~5 Myr 
old star. 

3. Previous modeling attempts to explain the structures 
in the disk fail, as noticed by the authors themselves, to 
reproduce them (Takeuchi & Artymowicz, 2001). The 
differential migration of dust grains in a gas disk with 
an outer sharp edge leads to the formation of a gap 
and an outer dust over-density of small particles but 
this approach, which apparently requires the presence 
of an unaccounted for sharp outer edge to the gas disk, 
results in a very axisymmetrical structure which does 
not match the more recent observations. 

2.2. Model and assumptions 

We consider in this paper the evolution of a coUisionless 
circumstellar disk of solid particles gravitationally per- 
turbed by an external bound perturber on an eccentric 
orbit. 

We assume that at least one of HD 141569 B and C and 

possibly both are boimd to HD 141569 A. Because the or- 
bital parameters of the two companions with respect to 
HD 141569 A are unknown and because of their close pro- 
jected positions we will adopt the simplest assumption 
that there is a single perturber assumed to be coplanar 
with the disk that can interact strongly with it. In princi- 
ple this could be either one of the companions or a compos- 
ite of them if they are bound together. We set the mass of 
the single perturber to the total mass of HD 141569 B and 
C assumed to be respectively M2 and M4 stars. But note 
self-similar properties of the simulations enable a scaling 
to different companion masses. This leads to a secondary 
to primary mass ratio of ^0.2. The pcriccntcr distance of 
the perturber is constant in the model and we vary the 
eccentricity of its orbit (Figure 2). Thus, unless explicitly 
stated, distances will be expressed in units of pericentre 
distance between the perturber and the primary. Four per- 
turber eccentricities are explored in this paper : e = 0.1, 
0.3, 0.5 and 0.7. The perturber starts the simulation at its 
pericentre position. 

The initial circumprimary disk consists of 9 x 10^ 
"massless" test particles distributed in a 2D axisymmetri- 
cal disk from r = l/6tor = 7/12 model units with a radial 
power-law surface density S(r) oc r". The disk is splited 
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Fig. 2. A basic sketch of the circumprimary disk at the begin- 
ning of a simulation and of the four perturber orbits considered 
in the paper. Distances arc expressed in model units (pericen- 
tre distance = 1). The orientation of the pericentre, along the x 
axis in the direction a; > 0, is fixed throughout the simulations. 

up into nine concentric annuli of 10^ particles each with ra- 
dius following the relation ri+2 = {2r'^^^ - r^"+2)V(«+2) 
ensuring the continuity of the surface density from one an- 
nulus to another. A constant surface density [a = 0) is as- 
sumed but different initial surface densities were explored 
by an a posteriori processing. This actiially does not im- 
pact the general behavior described in this paper. Particles 
arc initially on circular orbits around HD 141569 A. 

The orbits of the particles are numerically integrated 
using a fifth-order Cash-Karp Runge-Kutta method. 
Effects on the motion arising from radiation are neglected. 
The simulations start with a time-step of a tenth of the 
orbital period of the closest particle to HD 141569 A. Then 
it is adjusted within the code in order to ensure accuracy. 
The dynamical response of the disk is numerically followed 
over a total span of 169, 116, 70 and 32.5 perturber orbital 
periods for the 0.1, 0.3, 0.5 and 0.7 perturber eccentrici- 
ties respectively. Every fiftieth of the total span, positions 
and velocities of the particles are stored. As an example, a 
pericentre distance of 1200 AU and an eccentricity e=0.5 
correspond to a total span of 5 Myr with a storage process 
time-scale of 10^ years. 

3. Numerical results 

Figures 3a to 3d show snapshots of the temporal evolu- 
tion of an initially axisymmetrical circumprimary disk of 
test particles in an eccentric binary system with a mass 
ratio of 0.2 for different perturber eccentricities. Several 
features characterize the dynamical response of the disk 
to the perturber among which the clearing out of the outer 
disk and the formation of spiral structures are the most 
prominent ones. 

Particles initially at radii larger than r = 0.4 — 0.45 
model units reach highly eccentric orbits after about ten 
perturber periods resulting in a gravitational truncation 
of the disk consistent with theoretical calculations (e.g. 
Papaloizou & Pringle, 1977). The truncation occurs re- 
gardless of the assumed perturber eccentricity e (Figure 
4) but low perturber eccentricities increase the efficiency 
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Fig. 3. Evolution of an initially axisymmctric circumprimary disk within a binary stellar system for different perturber eccen- 
tricities (see section 2). The images show face-on views of the disk surface density in a logarithmic scale. The upper left image 
shows the common initial disk with 9 x 10^ particles spread over a radial range from r = 1/6 to r = 7/12. For each perturber 
eccentricity, relevant snapshots have been selected in order to highlight the dynamical response of the disk to the gravitational 
perturbation of the perturber. For each panel, the time is expressed in perturber orbital periods. The position of pericentre of 
the perturber is along the x axis in the right direction (see the sketch in Figure 2) . The position of the central star is indicated by 
the white cross. Brightness levels are the same from one image to another allowing direct comparison between all the snapshots. 



of particle ejection after a fixed number of orbital peri- 
ods since the pericentre distance of the pertnrber is con- 
stant in the model. For instance 90% of the particles ini- 
tially outside the critical radial distance of r = 0.45 model 
units have been placed on unbound orbits for e = 0.1 af- 
ter 10 perturber revolutions whereas 45% of them were 
still bound to the central star for e = 0.7 after the same 
number of orbits (Figure 4). Together with the low eccen- 
tricities induced by the perturber on the particles interior 
to r = 0.4 model units, this effect results in sharp outer 
disk edges for low e values. Conversely the disk truncation 
is less marked for the largest e values. This is contributed 
to by the increase of the mean eccentricity of the bound 
particles which have the ability to perform excursions sig- 
nificantly beyond the truncation distance (Figure 4b and 
c). Actually the disk truncation for large e values is not 



homogeneous in azimuth as can be noticed in Figure 3 for 
e = 0.5 and e = 0.7. This is manifest by the formation of 
a marked azimuthal asymmetry of the surface density for 
distances larger than r = 0.4 — 0.45 model units for the 
largest e values after similar numbers of perturber orbital 
periods. 

Inside the truncation distance, an over-density of par- 
ticles develops in conjunction with the truncation process 
which breaks the axisymmetry of the initial disk. This spa- 
tially coherent over-density appears close the outer edge 
of the truncated disk at a distance r ~ 0.35 model units. 
As shown below, the secular perturbation of the disk by 
the perturber is essentially responsible for the formation 
of this spiral-shaped structure which is well marked for 
e > 0.3. It progressively winds around the central star 
becoming radially thinner with time. The formation of 
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e = 0.1 6 = 0.3 e = 0.5 e = 0.7 




Fig. 4. Panels (a) and (b) : truncation distance versus time. Distances outside which the surface density is less than half 
(rcsp. a quarter of) the initial surface density are represented in panel (a) (resp. (b)). Time is expressed in number of perturber 
orbital periods. Distances are expressed in model units, i.e. normalized to the perturber pericentre distance. Panels (c) : 
mean eccentricities (lower panels) and fractions (upper panels) of bound particles versus time for the four different perturber 
eccentricities and four ranges of distances to the primary star (in units of perturber pericentre distance) . Time is expressed in 
number of perturber orbital periods but the total physical time is the same for each value of perturber eccentricity (e). The 
range of distances represented by the solid line on each panel brackets the truncation distance of the disk. 



the spiral-shaped over-density of particles leads in turn to 
the formation of lower surface density regions and ulti- 
mately results in the opening of a spiral gap. The gap 
always appears to be radially wider than the adjacent 
over-density spiral. The contrast between the high and 
low density regions begins to increase with time as the 
spiral propagates inwards. The particle eccentricity in- 
duced by the perturber depends on its eccentricity (Figure 
4). This prevents the spiral structure from having exact 
self-similar properties as the perturber eccentricity varies. 
Nevertheless we notice that the large scale features evolve 
in a similar way as a function of the perturber orbital pe- 
riod (Figures 3b to 3d). The latter characteristic is used in 
section 5 in order to derive orbital parameters for the per- 
turber matching the observations of the HD 141569 sys- 
tem. 

4. Spiral structure by secular perturbation of the 
inner disk 

We here show that the transient spiral structure in the 
inner disk seen in the simulations can be understood as 
being due to the action of the time averaged potential due 
to the perturber. This accounts for some of the self-similar 
scaling properties seen in the simulations. 

We adopt a simple model based on calculating the 
linear response of a cold collisionless disk to an orbiting 



companion in the continuum limit. In this scheme differ- 
ential precession would cause the induced spiral structure 
to wind up indefinitely generating arbitrarily small scales 
and large surface densities. In a more realistic situation, 
the fact that the disk is not a continuum and has a non 
zero velocity dispersion will cause the spiral structure to 
eventually wash out. To model such effects we introduce 
an ad hoc decay rate for the spiral form that increases as 
its radial scale decreases. Further we impose a cut off on 
the magnitude of the surface density calculated using the 
linear response. In this way the temporal morphology of 
the spiral form apparent in the linear response remains 
while unphysical effects arising through the simplicity of 
the model are suppressed. 

Wc suppose that (r, ip) define a cylindrical coordinate 
system based on the primary star of mass M* which exerts 
a potential ^ = —GMt,/r. We take the perturber with 
mass Mp to be in an eccentric orbit with semi-major axis A 
and eccentricity e. The apsidal line is taken to lie along the 
X axis {if = 0). The perturbing potential ^'p it produces 
can be expanded in terms of a Fourier cosine series in ip 
and time averaged (eg. Terquem & Papaloizou, 2002). It is 
necessary to retain only the first two terms in the Fourier 
series to lowest order in r so that we may write 

GMpr^ fl 3ercos^ \ 
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Fig. 5. Surface density plots taken at the number of orbital periods after initiation indicated. These plots are for a pertuber 
with e = 0.7 and Mp = 0.2M* . The analytical procedure described in section 4 was used. The pericentre distance was set to 1 

distance unit and the direction of pericentre is the same as in Figures 2 and 3. Disk radii between 1/6 and 0.47 arc considered 
which can thus be directly compared with Figure 3d. Such a comparison shows that similar morphology develops at similar 
times. 



The first axisymmetric term (m = 0) , ^po, contributes to 

the unperturbed motion in the disk. In particular it leads 
to a precession of orbits with small eccentricity at a rate 
u}p = fl — K, with O and k being the angular and epicyclic 
frequencies associated with a circular orbit containing the 
conserved orbital angular momentum respectively (see eg. 
Binney & Tremaine, 1987, for the basic definitions and 
disk dynamics required). For the time averaged potential 
\l/+\l/pO, assuming, as is the case, that uip <C ^, one obtains 



3GMp 



2A3(l-e2)3/2- 



(2) 



In general this is a function of the radius, ro, of the circular 
orbit. The presence of the second nonaxisymmetric (m = 
1) term, *pi in (1) perturbs these to become eccentric. 
One may write 



(3) 



where (^^o ) Cvo ) denote the components of the Lagrangian 
displacement vector in the radial and azimuthal directions 
respectively. Note too that (po = ^iro)t + /3o, where /3o is 
a fixed phase associated with and identifying a particular 
orbit. 

Assuming = ^d^o, with being a small disk eccen- 
tricity, in the linear approximation the governing equation 
is that of a forced harmonic oscilator taking the form 



dr r 



(4) 



r=ro,V=V>o 



Note that the time derivative is taken along a given unper- 
turbed orbit and so evaluation of the perturbing potential 
at y = (^0 introduces a time dependence on such an orbit. 

Equation (4) is easily solved for provided appro- 
priate boundary conditions are applied. Here we assume 
that prior to some time t = Q when the companion is 
introduced the disk is unperturbed. It is important to em- 
phasize that the formation of transient spiral structure re- 
quires such a 'sudden' introduction of the perturber. Here 
'sudden' means fast compared to the inverse precession 



frequency at the disk location of interest. A possibility 

would be a distant scattering into an orbit of high eccen- 
tricity near apocentre. A very 'slow' introduction of the 
perturber through build up of its mass on a long timescale 
would not lead to transient spiral structure. A solution of 
(4) corresponding to 'sudden' introduction in the limit of 
small ujp/Q. can be expressed in the form 



4A{1 - e2) 



(cos((y9o) - cos((yPo - Wp(ro)t)) . (5) 



The first thing we point out about (5) is that it has a 
self-similar scaling property for e ^ 1. Noting that in that 
limit A{1 — e^) ~ 2A{1 — e) being twice the pericentre 
distance and using equation (2), we see that equation(5) 
gives the same displacement for a given pericentre dis- 
tance and perturber mass ratio after a given number of 
orbital periods independent of A. This scaling, which is 
only appropriate for secular perturbations, occurs in our 
simulations. 

Note too that in the expression (5) the perturber 
mass only occurs as linear multiple in tUp as given by 
equation(2). Thus different perturber masses produce the 
same spiral planform, but on longer timescales for lower 
mass perturbers. We have also verified that the simula- 
tions show this. Adjusting either one of the companion 
masses or their combination can then be accommodated 
by scaling the time. 

Another aspect indicated above is that because Up de- 
pends on ro we have differential precession and an un- 
limited winding up and shortening of the scale of the 
disk perturbation. This occurs through an increase of the 
wavenumber k = {du!p/dro)t that is implied in the second 
term in (5) which is there because of the 'sudden' initial 
condition. Such large values of k can in theory lead to un- 
realistically large or even singular values for the surface 
density. 

In fact an arbitrary increase of k does not manifest it- 
self in the simulations or the physical problem of interest. 
It does so in the simplified analysis given above because 
the disk was assumed to be a continuum initially with ex- 
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Fig. 6. Left panel: dust surface density derived from HST/ACS images after deprojection, correction for stellar flux dilution 
and correction for anisotropic scattering (Figure 7d from Clampin et al., 2003). Middle and right panels: simulated disk surface 
density after ~9 perturber orbital periods. For illustrative purpose, an argument of pericentre of cj = 50° has been chosen. The 
true anomalies of HD 141569 B and C are indicated on the right panel assuming a disk inclination of 55° from pole-on and an 
observed disk position angle of -3°. This leads to deprojected distances for HD 141569 B and C of 1078.17 AU and 1291.29 AU 
respectively. 



act circular orbits only. The introduction of some velocity 
dispersion would result in a particle sampling a range of 
radii and thus a limit to the size of the value k and the 
disk surface density that may occur in any visible pertur- 
bation. Very short scales, if initiated, would be washed out 
and decay. Here we represent the decay of very short scale 
disturbances in a simple ad hoc manner. We simply reduce 
the amplitude of the cos((y9o — LUp{rQ)t) term in equation 
(5) by a factor that depends only on time by replacing 
it by cos(</?o — ^p{i'o)t)/ \/^ + 0.025k'^R^ /ro, where i? is a 
fiducial radius. This reduction causes the second term in 
(5) and hence the spiral structure to decay with time. It 
is clear that the first and remaining term at large times in 
(5) does not produce a spiral form. 

We now compare the form of the surface density per- 
turbation resulting from the above analysis with that 
found in the simulation. To do this we use mass conserva- 
tion expressed in Lagrangian form, namely 



where J = \d{r,ip)/d{ro,'Po)\ is the magnitude Jacobian 
determinant of the transformation {r,ip) (ro,</Jo)- Here 
S is the surface density of the disk under perturbation, 
while So is the surface density for the unperturbed disk. 
We adopt the relation between s-nd that ap- 
plies to Keplerian orbits with small eccentricity, namely 
d^,^„/dipQ = — 2^ro- To prevent unphysically large surface 
densities we employ a cut off through providing a floor for 
J which limits the surface density at 2So. 

We plot the surface density obtained from the above 
analysis in Figure 5 for a binary with e = 0.7, Mp = 0.2Mh< 
and pericentre distance 1 model units. The disk is consid- 
ered for 1/6 < ro < 0.47 model units and that corresponds 



to the section of the disk illustrated in Figure 3. The num- 
ber associated with each plot is the time in orbital periods. 
Comparison with Figure 3d for the simulations indicates a 
similar morphology regarding the spiral pattern at corre- 
sponding times. Note that we have verified that this mor- 
phology is insensitive to the mode of the cut off at high 
surface densities or of large radial gradients. Furthermore 
because the morphology is generated by the time averaged 
secular potential, it does not depend on the precise loca- 
tion of the perturber on its orbit at any time but only on 
the direction to pericentre, here being along the x axis. 

5. Application to the HD 141569 A disk 

We discuss in this section the morphological similarities 
between the disk features described in the two previous 
sections and those observed in the HD 141569 A disk. We 
base our semi-qualitative comparison between simulations 
and observations on the dust surface density maps ob- 
tained by Clampin et al. (2003, hereafter COS). These 
maps have been derived from HST/ACS deprojected im- 
ages of the disk seen in scattered light in B, V and I bands. 
They take into account the stellar flux dilution effect as- 
suming the disk is optically thin in all directions at visible 
wavelengths (Figure 7c from C03). They take also into ac- 
count the anisotropic scattering properties of the grains 
assuming forward scattering with an asymmetry factor g 
in the range 0.25-0.35 consistent with the upper limit as- 
sessed by Mouillet et al. (2001) (Figure 7d from C03). 
Given the relative uncertainties in grain properties in the 
HD 141569 A disk (see nevertheless section 6), and thus 
the uncertainties on grain scattering properties, we may 
safely assume in the following that Figures 7c and 7d in 
COS bracket a representation of the true surface density of 
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the dust disk by noting that they respectively correspond 
to the isotropic and maximum anisotropic assumptions for 
scattering. 

We interpret the following observational evidences : 

— the lack of obvious link between the resolved inner and 
outer annulus in the HD 141569 A system, 

— the low surface density of the inner ring compared to 
the outer one, 

— the azimuthal asymmetry of the outer ring, 

as three constraints on the number of perturber orbital 

periods the system has evolved for and on the eccentricity 
of the perturber's orbit. According to these criteria, well 
wound structures are ruled out such as those formed af- 
ter '~20 perturber orbits. The "isotropic" surface density 
is better recovered after ~15 perturber orbits as against 
^10 perturber orbits for the "anisotropic" case (Figures 
3b to 3d) . The measured azimuthal asymmetry within the 
structure is reproduced for e > 0.3 (Figure 7). For e = 0.1, 
the disk is too poorly contrasted both azimuthally (Figure 
7) and radially (Figure 3). 

We then conclude that the surface densities of 
HD 141569 A inferred from the HST/ACS observations 
show large scale features consistent with a circumprimary 
disk within an eccentric binary system provided that the 
perturber eccentricity is at least larger than 0.3 (most 
probably 0.5) and second that the system has been per- 
turbed for a total duration of ~ 10-15 orbital periods 
(Figure 6). Under such conditions the spiral-shaped struc- 
ture newly-formed close to the truncation distance quali- 
tatively matches the predominant structure observed at a 
distance of 325 AU (also highlighted in red color on Figure 
8a from C03). The position of the over-density which ap- 
pears close to r = 0.35 model units after ^10 15 per- 
turber orbits (sections 3 and 4) matches the observations 
if the pericentre distance is set to ^930 AU. An argu- 
ment of pericentre of ^50° (anti-clockwisc direction on 
the Figures) with an uncertainty of at least 10° qualita- 
tively matches the observations. Such orbital parameters 
for the perturber are consistent with current deprojected 
positions of the two companions assumed to lie in the disk 
plane and assimilated into a single perturber in our ap- 
proach (Figure 6). 

In a direction opposite to that of the over-density of 
particles, the simulated disk shows a broad faint exten- 
sion. A first consequence is that the outer edge of the 
disk looks shifted in the direction of this extension while 
the central disk hole (artificially introduced in the sim- 
ulations) remains centered onto the star. Therefore, the 
formation of a spiral structure in the outer regions of the 
disk as well as the shift of the center of the outer disk edge 
provide explanations of the shifts reported by several au- 
thors between the inner and outer "rings" (Mouillet et 
al., 2001; Boccaletti et al, 2003; Clampin et al., 2003). 
A second consequence is that the broad extension oppo- 
site to the over-density may account for the "diffuse emis- 
sion" reported by Mouillet et al. (2001) in that direction 
(see Figure 1). An alternative explanation for the "diffuse 
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Fig. 7. Azimuthal asymmetry versus time averaged in a radial 
region from r — 0.33 to r = 0.37 and for different perturber 
eccentricities (e). The azimuthal asymmetry is here defined as 
the ratio of maximal to minimal total numbers of particles in 
30° angular sectors. It has boon derived from the simulations 
presented in Figure 3. The numerical noise leads to a typical la 
uncertainty of 0.1. The radial range considered here includes 
the over-density of particles appearing after ~10 perturber or- 
bital periods and which compares well with the observations 
(see section 5). The brightness asymmetry measured in the vis- 
ible by Clampin et al. (2003) is indicated by the light grey area. 
A perturber eccentricity e > 0.3 and 10-15 perturber orbital 
periods arc both required in order to reproduce the observed 
brightness asymmetry. 



emission" could be a local enhancement of the collision 
frequency at the location of the over-density resulting in 
the production of a larger number of small particles with 
large eccentricities. 

Matching the estimated star age of 5 ± 3Myr 
(Weinberger et al.. 2000; Merfn et al., 2003) critically re- 
duces the range of companion eccentricities : a 2 Myr age 
implies e = 0.78±0.03 while 8 Myr leads to e = 0.91±0.01. 
We made use of the rough self-similar properties applying 
inside the truncation distance as noted in sections 3 and 
4 in order to derive these values. Lower eccentricities are 
theoretically possible but less probable in the framework 
of a gas-free disk of solid material unless the perturbation 
is very recent (see the discussion in section 4). 

6. Minimum grain size 

Clampin et al. (2003) also report the redness of the disk 
in the visible. The disk scatters 25 ± 2% more stellar flux 
in the I band than in the B band and 10 ± 0.7% more 
in the I band than in the V band. Interestingly no color 
gradient within the disk is measured which indicates that 
similar scatterers dominate the scattered light images at 
every distance from the star. One can anticipate that the 
grains responsible for the color effect are close in size to 
the wavelengths discussed here. In this section we further 
constrain the size distribution in the disk and discuss some 
implications. 
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Table 1. Constraints on the minimum grain size Umin in the 
disk based on the interpretation of the color indexes mea- 
sured in the visible (HST/ACS, Clampin et al., 2003). Compact 
grains with a differential size distribution proportional to a"" 
arc assumed. Sizes are expressed in fim. They are computed 
with Eqs. A. 8, A. 11 (with a = 3) and A.l assuming k = 3.5 and 
complex dielectric permitivities of e = 2.8 + iO.l (or n ~ 1.67) 
for amorphous silicate and e = 5.1 + iO.6 (or n ~ 2.26) for 
graphite (Laor & Draine, 1993). "B-I" and "V-I" refer to the 
bands used to measure the color effects, oo is computed in the 
I band. 





lower limit 


upper limit 




on ttmin 


(amin) 


on ttmin (10 ao) 




B-I 


V-I 




Silicate : 


0.077 


0.095 


3.1 


Graphite : 


0.041 


0.051 


1.6 



than the blow-out size limit (so-called /3-meteoroids in the 
Solar System). A similar result was derived from the in- 
terpretation of a marginal color eflfect in the near-infrared 
(Boccaletti et al., 2003). 

According to the results summarized in Table 1, the 
size of the smallest grains in the HD 141569 A compares to 
the largest interstellar grains (e.g. Li & Greenberg, 1997; 
Weingartner & Draine, 2001; Clayton et al., 2003), to the 
monomers of comctary dust particles (e.g. Greenberg & 
Hage, 1990; Kimura, Kolokolova, & Mann, 2003) or to 
the constituent monomers of interplanetary dust parti- 
cles, thought to be of comctary origin, and which exhibit 
porous aggregate structures (e.g. Brownlee, 1985). The in- 
terpretation of the scattered light observations leads to 
minimum grain sizes consistent with that proposed by Li 
& Lunine (2003) from their analysis of the thermal emis- 
sion of the disk. They show that the spectral energy distri- 
bution of the disk is well reproduced if 0.1 < amin < 1 M™- 



6.1. Interpretation of the disk color in the visible 

For simplicity we assume that the differential dust size dis- 
tribution dn(a) follows a single power law; dn{a) cx a^'^da 
between a minimum size ftmin and a maximum size Omax- 
We in addition restrict the discussion to compact grains""^. 
The HD 141569 A dust disk shares some similarities with 
typical debris disks: it is optically thin and grains are 
short-lived (Boccaletti et al., 2003; Li & Lunine, 2003). 
A reservoir of mass, presumably an unseen population of 
colliding large bodies such as planetesimals, feeds the disk 
in observable dust grains. The dust size distribution there- 
fore extends to large particles that behave as grey scatter- 
ers in the visible and near-IR. The color effect measured 
in the visible is very likely produced by the smallest par- 
ticles of the size distribution which implies k > 3 (mean 
scattering cross section dominated by the smallest grains). 

In the small range of wavelengths considered here, a 
noticeable color index is inconsistent with dominant scat- 
terers in the regime of geometric optics. Grains larger than 
about ten times the size ao at which the scattering effi- 
ciency Qsca reaches its maximum induce fluctuations of 
<3sca of the same order of or smaller than the uncertainty 
on the color index reported by C03. This remark leads 
to very conservative upper limits on amin- On the other 
hand a strict lower limit on amin consistent with the obser- 
vations can be derived by considering the grain scattering 
behavior in the Rayleigh regime. This approach is detailed 
in Appendix A (Eq. A. 11). Table 1 gives the results for 
spherical silicate and graphite compact grains. 

The upper limit on the minimum grain size in the disk 
must be compared to the blow-out size Op,. in the radiative 
environment of HD 141569 A. Boccaletti et al. (2003) show 
that compact spherical grains smaller than Upr ~ 6 fim are 
blown away by radiation pressure. We conclude from our 
interpretation of the redness of the disk reported by C03 
that the disk contains a large fraction of grains smaller 

^ The impact of the porosity is discussed in section A. 3 



6.2. Consequences 

The presence of a large amount of small grains, as small 
as a tenth of the blow-out size limit, points out on the 
need for an efficient mechanism able to continuously re- 
plenish the dust disk in fresh particles, theoretically ex- 
pelled from the disk after a few orbital periods, and re- 
sponsible for the observed color effect. Mean collision 
time-scales two or three orders of magnitude smaller than 
the star age have been estimated from the vertical opti- 
cal thickness of the HD 141569 A disk (Boccaletti et al., 
2003). Therefore coUisional evolution can happen. Given 
the grain sizes and the large distances from the star con- 
sidered here, Poynting-Robertson drag is not efficient com- 
pared to collisions. The large number of /3-meteoroids in 
the disk, large enough to produce a noticeable color effect, 
probably results from the large collisional activity of the 
disk compared to more evolved debris disks around older 
stars. According to Krivov, Mann, & Krivova (2000) and 
Augereau et al. (2001), the outer disk of /JPictoris could 
also hold a population of /3-meteoroids. But the shape 
of the size distribution departs from a single power law 
in Krivov. Mann, & Krivova (2000) 's model^ and a ma- 
jor consequence is that the cross-section is dominated by 
bound particles that have a grey scattering behavior con- 
sistent with the neutral color of the outer /?Pictoris disk 
in the visible (e.g. Lecavelier Des Etangs et al., 1993). 

Alternatively, we can further consider the optical be- 
havior of the grains. The scattering properties of an aggre- 
gate might be similar to that of the constituent particles if 
the aggregate is porous enough. Could the dominant scat- 
terers responsible for the images of the HD 141569 A disk 
be monomers packed into porous bound grains? We let P 
the porosity of an aggregate and ap^{P) the correspond- 

^ Using a different approach, Thebault, Augereau, & Beust 
(2003) also show that dust production by collisions results in 
size distributions that significantly depart from single power 
laws 
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ing blow-out size limit. For large P values, the number of 
monomers required to form a grain of size apr(f ) can be 
approximated by 



iV ~ (1 - P) 



pr(P) 



(1-P)- 



apr(P = 0) 



where is the mean size of the monomers. For 
HD 141569 A, the ratio Opr(P = 0)/am is of the order of 10. 
If the scattered light images are indeed tracing monomers 
in very porous bound aggregates, then the smallest grains 
in the HD 141569 A system are made of at least a few thou- 
sands of monomers similar in size to the constituents of 
the interplanetary dust particles. The major problem with 
this alternative scenario is that collisional time-scales are 
very likely short enough to efficiently split up such large 
aggregates into their constituent monomers producing an 
unavoidable population of particles small compared to the 
blow-out limit. 

7. Conclusion 

We have explored in this paper the possible impact of the 
two M companions of HD 141569 A on the dynamics of its 

circumprimary dust disk. Provided that at least one of the 
two companions is bound to HD 141569 A and on an orbit 
of high eccentricity, the tightly wound and asymmetric spi- 
ral structure at ~ 325 AU can be reproduced. By contrast, 
the disk interior to the spiral structure looks depleted as 
observed. The inner ring which has a low surface density 
compared to the outer structure is not reproduced in this 
approach. From the redness of the disk in the visible we 
deduce that the disk contains a large fraction of short- 
lived grains small compared to the blow-out size limit and 
implying high collisional activity. 

How the contrasting structures of the HD 141569 A 
disk are affected by collisional activity and the effects of 
radiation pressure are key issues that will need to be ex- 
plored. The dynamics of the disk is made even more com- 
plex by the fact that gas drag may also play a role. The 
system indeed contains a remnant amount of gas that was 
detected for the first time by Zuckerman et al. (1995). 
Whether the gas disk is extended and massive enough to 
impact the dynamics and the shape of the dust disk es- 
pecially at very large distances (around 300-350 AU) is an 
issue that will be addressed in a near future thanks to re- 
cently resolved images of the CO disk (Augereau et al., 
2003). 

Our model does not provide an explanation of the ob- 
served depletion of solid material inside 150 AU. An addi- 
tional process is required to clear out the inner disk and to 
produce the sharp edge observed around 150 AU. Together 
with the surprising detection of H^ in the HD 141569 A 
disk by Brittain & Rettig (2002), previously observed in 
atmospheres of the giant Solar planets, these observa- 
tions are suggestive clues for speculating on the presence 
of planetary companions in the inner disk. The coupling 
of putative newly-formed massive companions in the in- 
ner disk with the gas disk resolved by Augereau et al. 



(2003) is an exciting prospect. The disk-planet coupling 
has been theoretically addressed by several authors and 
the orbital parameters of some of the detected extra-solar 
planets are understood as the final result of this interac- 
tion. The HD 141569 A system is young enough (only a 
few Myr) with a still reasonably large amount of remnant 
gas so that disk-planet coupling might still act. Future ob- 
servations of the HD 141569 A system should judiciously 
focus on the inner regions to better assess how empty and 
perhaps structured the inner gas and dust disk is, and 
on opportunities for detecting possible sub-stellar objects 
inside 150 AU. 

Acknowledgements. Wc thank A.M. Lagrange and D. MouiUct 
for helpful discussions on the interpretation of the high res- 
olution images of the HD 141569 A disk. We also thank M. 
Wyatt for carefully refereeing the paper, P. Thebault for use- 
ful comments and the ACS Science Team for kindly provid- 
ing HST/ACS images. J.C.A. is deeply grateful to E.F. van 
Dishoeck for her permanent support in the course of this work 
at Leiden Observatory and to C. Eiroa for the discussions on 
the distance of HD 141569 A. J.C.A. was supported by a fel- 
lowship from the European Research Training Network "The 
Origin of Planetary Systems" (PLANETS, contract number 
HPRN-CT-2002-00308) at Leiden Observatory. 

Appendix A: Interpretation of a scattering color 
effect in terms of minimum grain size 

The Mie scattering efficiency Qsca of spherical and ho- 
mogeneous dust particles only depends on the real part 
n\ = Re{m\) of the complex refraction index of the grain 
fnx = \/eA and on the so-called size-parameter x = lnajX 
where a is the grain radius and A the wavelength. Q^ca. fol- 
lows three basic well-known regimes (Bohren & Huffman, 
1983) : 

1. for small x values (a <C A, Rayleigh regime), Qsca oc a;^. 

2. for large x values (a ^ A) Qsca ^ constant (grey scat- 
tering). 

3. otherwise incident rays and rays passing through the 
grain produce interferences resulting in extinction 
maxima (then scattering maxima) peaking at Xi = 
{2i + l)7r/2(nA — 1) where « is a positive integer or 
zero. For a given wavelength, the first peak appears 
for grains of size 



ao 



X 



4(nA-l) 
We note 

xo = 2nao/X 



• - 1) 



(A.1) 



(A.2) 



We propose an approach allowing the interpretation of a 
measured scattering color effect in the context of an opti- 
cally thin circumstellar disk. The color effect is expressed 
by 



C 



$2/^2 
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where <5?i and $2 are the measured circumstellar flux at 
wavelengths Ai and A2 > Ai; and $2 are the stellar 
flux at Ai and A2. Wc assume in the following that the 
grains have similar complex optical indices in the range 
of wavelengths considered. We note n = n\ and m = m\. 
Anisotropic scattering properties arc then similar and the 
scattering efficiency Qsca only depends on the size param- 
eter X = 27ra/A for spheres of radius a. In the optically 
thin regime the color effect can be then written 



c = 



where (cri) and ((72) are the scattering cross-sections aver- 
aged over the grain size distribution dn{a) cx a^'^da. a,nin 
is the minimum grain size and flmax is the maximum grain 
size. For the sake of simplicity we assume 3 < k < 7 . We 
also assume that amax ^ ao at both wavelengths Ai and 
A2. Thus, the color effect does not depcind on Omax and 
this allows us to replace Omax by +00 in the following. 
The color effect C can be rewritten 



C = 



f 



+ 00 



(x)A?-'^a;2-«da; 



where 



Xi 



and X2 = 



27ran 



Ai A2 
Since Xi > X2 we rewrite C in the form 



K — 3 



rx 
J X2 



(A.3) 



A.l. Upper limit on Omin in the Rayleigli regime 

We now consider minimum grain sizes a^i^ such that the 

scattering efficiencies Qsca.{xi) and Qsca.{x2) fall into the 
Rayleigh regime. The basic idea is to compute an upper 
limit of the minimum grain size in the circumstellar disk 
consistent with the observed color effect in the regime of 
grains small compared to the wavelength. This gives 



/" 



i(a;2) 



1 



with 

Qsca{xo) < Qsca{X2) 



(A.7) 



By combining equations A.3 and A. 4 and inequalities A. 5 
to A.7 we finally obtain 



A I (A2/A1) 



7-K 



1 \ 



l/(«-7) 



X2 < Xq 



1-C(A2/Ai 



1 + 



7-K 



This leads to amin < aZ,;„ where 



A2 



A I (A2/Aif--1 \ 



l/(«-7) 



4(n- 1) 



1-C(A2/Ai) 



V 



1 + 



7-K 
K-3 



. (A.8) 



Minimum grain sizes amin obeying to the Rayleigh ap- 
proximation and consistent with the color measurement 
C must then be smaller than a^^^. We recall that the for- 
mer equation does ensure that such grains do exist in the 
disk. Larger grains, also consistent with the color index 
but that do not follow the Rayleigh approximation are 
not excluded. 

A. 2. Lower limit on Omin the Rayleigh regime 

A lower limit on the minimum grain size in the Rayleigh 
regime and consistent with the observations can similarly 
be derived provided that a real a fulfilling the following 
criteria can be found : 

1. 1 < a < xo/x2, 

2. Qsci,ixo/a) = Qsca(a;2)/a;2 x (xo/a)'', i.e. Qsca(a;) 
obeys the Rayleigh approximation for x from X2 to 
xo/a, 

3. Qsca(a^o/Q:) < 1, i.e. the scattering efficiency for x = 
xo/a is smaller than the value of Qsca{x) in the geo- 
metric optics regime. 



-tA-4) 



We compute the integral between X2 and -|-oo in equation 
A.3 in two steps. Between X2 and xq (given by equation 
A. 2) the Rayleigh regime is accurately obeyed except for 
X values close to Xq where the Rayleigh law over-estimates 
Qsca- Therefore 



Then 

Xo/a 



I 



Qsc3.{x)x'^ ''dx = 



1X2 

and 



Q5ca(x2) 

x5{7-K) 



(^)'-"-xrtA.9) 



f 



Qsca(a:2) 

x4(7-k) 



(xq a;2 



(A.5) 



The scattering efficiency reaches its maximum value for 
X = Xq. Assuming that Qsca is a constant for x larger 
than xq leads to an upper limit of the integral between xo 
and -|-oo 



'xo/a 

or 

xa/a 



/ ^sca(^)'^ dx ^ / Qsca ( 
Jxq/ol J xo/a ^ 



"-dx 



Q.Mx^-^dx>9if^(^)'-' 

X2{K — 6) \ a / 



(A.IO) 



/■+00 



Inequality A.IO together with equations A.3, A. 9, A. 2 and 
A.8 give 



Qsca.{x)x^ "^dx < Qsca{Xo)- 



(A.6) 



> a^i„ where a^-^ 



Ja . 



(A.ll) 



12 



Augereau & Papaloizou: Structuring the HD 141569 A circumstellax dust disk 



We add a few useful remarks about the criteria defin- 
ing a. The first criterion must be used to check the self- 
consistency of the result. The second criterion can consider 
the fact that size parameters a few times smaller than xq 
already ensures that the scattering efficiency follows the 
Raylcigh law with a good approximation. In order to ful- 
fill the third criterion one should use the precise Rayleigh 
approximation (e.g. Bohren & Huffman, 1983) 

2 



1 



e + 2 



Solving the previous equation for x ~ x^/a and Qsca.{x) = 
1 gives an upper limit an a. If grains are very porous and 
made of inclusions small compared to the wavelength, then 
e should be replaced by eoff where eoff is given by equation 
A. 12 and xq is given by equation A. 14. 

A.3. Grain porosity 

Equations A. 8 and A. 11 which depend on ao (or xq) have 

been obtained by assuming that grains arc spherical and 
compact. In the following, we assess the impact of the 
grain porosity P on ao (or xq). A way to go towards the 
solution is to use an effective medium theory as described 
in details in Bohren & Huffman (1983). Here we consider 
the Max;well-Garnett theory. In that framework, a porous 
grain made of a single material is likened to a two compo- 
nent grain made of a matrix of vacuum and inclusions of 
solid material. It is important to note that this theory is 
only valid for inclusions (considered as inhomogencities 
into the matrix) small compared to the wavelength. If 
these conditions are fulfilled then the effective dielectric 
permettivity of the porous grain is then given by 

3(l-P)(£-l)/(6 + 2) 



CeflF = 1 



(A.12) 



-(l-P)(e-l)/(e + 2) 

where e = er + iei is the complex dielectric permettivity 
of the solid material (inclusions in the Maxwell-Garnett 
approach). 

The real part of Ces can be approximated by n\ ~ 
^/Re{e^ provided that /m(ecir)^ Re{eesY. If the 
latter condition is fulfilled, then the Taylor expansion of 
{n\ — 1)/(1 — P) with respect to the porosity P gives 

^^ = (y^-l)+p,P + 0(p2) 

with < P < 1 and 



Pi 



Se,- - 6^67 + - eg 
6^ 



Therefore, the grain size for which Qsca reaches a first 
maximum is 



A 



ao = 



A(l-P)-i 



4(nA-l) 4(V^-l + piP) 
The corresponding size parameter is 
7r(l-P)-i 



Xq 



2(^-l + piP) 



(A.13) 



(A.14) 
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